National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


rrggrgiiS! 


*(  V 

I 


iiU 


u^Jb 


The  Rand  Paper  Series 

Papers  are  issued  by  The  Rand  Corporation  as  a service  to  its  professional  staff, 
Theii  purpose  is  to  facilitate  the  exchange  of  ideas  among  those  «vho  share  the 
author's  research  interests;  Papers  are  not  reports  prepared  in  fulfillment  of 
Rand's  contracts  or  grants.  Views  expressed  in  a Paper  are  the  author's  own,  and 
are  not  necessarily  shared  by  Rand  or  its  research  sponsors. 

The  Rand  Corporation 
Santa  Monica,  California  9040f> 


Far  Infrared  and  Submillimeter 


Wave  Attenuation  by  Clouds  and  Rain 
ABSTRACT 


Newly  determined  optical  constants  for  water  at  far  infrared  and 
subaillimeter  wavelengths,  as  revealed  by  a recent  survey, 'are  used 
to  estimate  water  cloud  and  rain  attenuation  over  the  wavelength 
range  between  X12pm  and  X2cm.  For  this  purpose  new  analytic  drop- 
?ize  distribution  models  simulating  fog,  nimbostratus  cloud,  and  rain 
corresponding  to  rainfall  rates  of  10  and  50  mm  hr  are  set  up.  The 
corresponding  volume  extinction  and  absorption  coefficients  are  computed 
according  to  polydi sperse  Mia  scattering  theory  at  specific  wavelengths 
and  presented  in  tables  and  graphically  in  plots  for  purposes  of  inter- 
polation. / 

It  is  found  that  cloud  extinction  may  exceed  50  nepers  per  kilometer 
in  the  A<100am  region  whereas  for  wavelengths  longer  than  X200pm,  under 
near  saturated  conditions,  water  vapor  absorption  should  be  the  dominant 
attenuator.  The  greatest  attenuation  by  heavy  rain  may  be  expected 
around  X5mm  with  a value  of  about  5 nepers  per  kilometer.  The  results 
also  suggest  that,  in  the  presence  of  non-precipitating  water  clouds 
or  fog  there  may  be  a relative  transmission  "window"  centered  around  XI. 3 
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ADDENDUM 


P-5419  FAR  INFRARED  AND  SUBMILLIMETER  WAVE  ATTENUATION  BY  CLOUDS  AND  RAIN, 
by  D.  Delrmendj ian,  April  1975. 


Please  add  the  following  as  a footnote  referring  to  the  word  "(neper)" 
appearing  on  line  5,  p.  15. 


As  used  here  and  elsewhere  in  this  paper,  the  term  "neper"  refers 
to  units  of  the  exponential  attenuation  of  power  or  energy  flux, 
rather  than  of  amplitude  or  voltage  as  in  accoustical  and  electrical 
engineering  usage. 
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1.  Introduction 

In  earlier  work  we  introduced  a set  of  continuous  size  distribution 
models  to  simulate  natural  atmospheric  particulates  such  as  aerosols  or 
haze,  c3oud  droplets,  and  raindrops  to  deduce  their  single  scattering 
and  absorption  properties  (Deirmendj ian,  1963a, 1964,  1965,  1969).  These 
models  have  been  found  useful  in  various  application  by  other  authors. 

The  work  covered  the  visible  and  infrared  spectrum  between  the  wave- 
lengths of  X0.45  and  X16.0  urn,  and  the  microwave  spectrum  between  XI  mm 
and  X8  cm  where  reasonably  accurate  data  on  the  optical  constants  of 
water — the  main  constituent  of  atmospheric  particles — were  available. 

As  is  well  known,  the  extinction  and  scattering  properties  of  individual 
particles  (van  de  Hulst,  1957)  as  well  as  of  polydispersions  (Deirmendj ian 

* 

1909)  depend  critically  on  the.  optical  constants,  i.e.,  the  real  and 

imaginary  parts  of  the  refractive  index,  of  the  particle  substance  with 
respect  to  the  ambient  medium. 

Two  recent  developments  make  it  poss.'ble  and  desirable  now  to  consider 
the  far  infrared  and  submillimeter  region  oetween  X10  pm  and  X1000  um: 

The  advent  of  several  types  of  new  maser  sources  of  monochromatic  radiation 
within  the  region  (Fox,  1971;  Button,  1974)  and  the  publication  of  new 
experimental  results  on  the  far  infrared  optical  constants  of  water  as 
reve.aled  by  a recent  survey  (Deirmendj  lan,  1974).  There  is  no  doubt  that 
these  new  sources  will  find  many  applications,  some  of  which  will  be 
related  to  the  terrestial  atmosphere  either  as  a subject  of  direct  study 
or  as  a disturbance  interposed  between  source  and  receiver.  Electro- 
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magnetic  wave  scattering,  is  connected  to  both  aspects:  as  a diagnostic 

of  particulates  through  the  study  of  the  directional  variation  of  scattered 
intensity  and  its  polarization,  and  an  an  attenuating  disturbance  through 
the  absorption  and  scattering  cross-section  of  the  particles.  As  is 
well  known,  the  most  efficient  scattering  is  produced  by  a particle  whose 
linear  dimensions  are  of  the  order  of  the  wavelength  of  the  incident 
radiation.  Cloud  and  rain  particles  have  typical  diameters  of  tens  and 
hundreds  of  micrometers  which  are  just  right  for  maximum  interaction  with 
far  infrared  (15  <X<  100  pm)  and  Bubmillimeter  (100<-\£l000  nm)  waves. 

The  present  paper  covers  mainly  the  attenuation  properties  of  these 
bydrometeors  leaving  the  directional  scattering  and  polarization  properties 
for  a projected  separate  study.  In  particular,  we  shall  first  introduce 
and  justify  some  new  cloud  and  rain  drop  size  distribution  models  developed 
for  use  in  this  study.  We  the.i  present  the  results  of  attenuation 
calculations  for  these  models  for  selected  wavelengths  and  optical  con- 
stants using  polydisperse  Mie  scattering  theory  and  computer  programs 
developed  earlier  (Deirmendj ian , 1969).  Finally  we  analyze  and  compare 
our  estimates  with  other  work  and  with  observation,  where  possible,  and 
discuss  seme  implications  in  atmospheric  applications. 

2.  The  optical  constants  of  liquid  water 

It  is  well  known  that  the  behavior  of  the  optical  constants  or  real 
and  imaginary  parts  v,  and  k,  respectively,  of  the  index  of  refraction 
u ; v * Ik,  for  liquid  water  is  rather  complex.  Particularly  in  the  middle 
and  far  infrared  region,  both  parts  are  highly  frequency-dependent,  and  vary 
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considerably  witi  temperature  in  a nonlinear,  nonmonotonic  way. 

Numerous  experimental  and  theoretical  techniques  and  models  are  used 

in  evaluating  these  so-called  constants.  Here  we  shall  merely  mention 

some  of  the  original  published  sources  we  have  used  in  the  present  work 

referring  the  reader  to  a longer  survey  for  more  details  (Deirmendjian,  1974). 

The  complex  index  of  refraction  m (with  respect  to  free  space)  is 

usually  defined  either  in  terms  of  Maxwellian  electromagnetic  theory  as 

the  square  root  of  a complex  dielectric  constant  or  a generalized  Snell's 

law,  provided  the  angle  of  refraction  in  a lossy  medium  is  also  defined 

as  a complex  number.  In  either  case,  the  real  part,  v may  be  related  to 

changes  in  the  wave's  phase  angle,  speed,  and  direction  of  propagation 

upon  entering  the  medium  from  free  space;  and  the  imaginary  part,  k, 

to  the  bulk  absorption  cross-section,  y,  , per  unit  volume  within  the 

b 

medium,  through  the  relation 


(1) 


is  nondimens tonal . For  weakly  abso-bing  or  "transparent"  media  the  real 
part  may  be  obtained  experimentally  by  an  appropriate  application  of 
Fresnel's  generalized  formulas  for  reflection  and  refraction.  This  method, 
of  course,  is  not  practical  for  water  in  the  far  infrared  because  of  its 
quasi  metallic  behavior  in  its  optical  properties,  and  other  experimental 
and  analytical  tools  must  be  devised. 

In  principle  the  imaginary  part,  k,  may  be  obtained  experimentally 
for  example  by  use  of  Beer's  lav;  by  measuring  the  absorption  through  a 
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variable  thickness  of  the  material.  In  practice,  however,  because  of  the 
very  high  infrared  absorptivity  of  water  one  has  to  use  a very  small  thick- 
ness— of  the  order  of  a few  wavelengths — of  the  sample  which  presents 
considerable  experimental  and  analytical  difficulties. 

The  optical  behavior  ot  liquid  water  varies  greatly  throughout  the 
spectrum,  from  that  of  a transparent  dielectric  in  the  visible  to  a 
metallic  type  in  the  far  infrared  region.  Some  of  the  difficulties  in 
measuring  the:  far  infrared  optical  constants  have  been  overcome  by  a 
combination  of  experimental  data  and  theoretical  models.  A brief 
description  of  these  methods  is  given  in  the  above  mentioned  survey 
(Del rmendj ian,  1974)  together  with  an  intercomparison  and  critique  of 
their  reliability.  Here  we  shall  merely  summarize  the  results  of  the 
survey  and  list  the  values  actually  used  in  our  calculations. 

Our  survey  indicated  that  the  optical  constants  of  water  in  the 
middle  and  far  infrared  below  about  A40  pm  are  fairly  well  determined. 

At  longer  wavelengths  there  are  uncertainties,  particularly  in  the  re- 
fractive index,  and  in  the  temperature  dependence  of  both  the  index  and 
the  absorption  coefficient.  It  is  possible  to  propose  a set  of  most 
probable  values,  based  on  the  latest  measurements  at  room  temperature, 
which  are  adequate  for  purposes  of  scattering  calculations.  The  results 
of  our  survey  are  conveniently  summarized  by  the  points  and  curves  of  Fig.  1 
plotted  on  logarithmic  scales  for  both  the  wavelength  and  the  optical  con- 
stants. The  ordinate  serves  a double  purpose:  the  left-hand  scale  is 

marked  in  units  of  the  imaginary  part  x , of  the  complex  index;  the  right- 


hand  scale  is  marked  in  units  of  (v  - 1)  for  convenience  of  representation. 
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The  continuous  curves  are  plots  of  the  smoothed  values  proposed  and  tabulated 
by  Hale  and  Querry  (1973)  and  based  mostly  on  work  by  a Kansas  State  Uni- 
versity group  initiated  by  Dudley  Williams  and  his  collaborators  (Draegert 
et  al,  1966).  In  Fig.  1,  the  solid  curves  for  water  ending  at  X50  pm  show 
their  "best"  values  and  should  be  regarded  as  fairly  illustrative — but  not 
definitive — for  the  region.  At  longer  wavelengths,  their  tabulations,  in- 
dicated by  the  dotted  curves  ending  at  X2Q0  pm,  may  be  less  tellable  and 
subject  to  future  correction. 

The  open  dots  in  Fig.  1 represent  the  optical  constants  adapted  from 
the  data  of  Davies  et  al.  (1970)  in  the  region  100  < X < 1600  pm.  In 
general  we  believe  that  these  values  represent  the  most  reliable  data  within 
their  range  of  coverage.  In  the  millimeter  region,  the  Davies  et  al.  (1970) 
experimental  values  at  XI  mm  compare  very  well  with  those  given  by  the 
Debye  formula  as  adapted  for  example  by  Ray  (1972)  for  water  at  30°C. 

For  completeness  and  comparison  we  have  included  in  Fig.  1 the  optical 
constants  for  ice  according  to  some  recent  data  (Deirmendj lan,  1979), 
in  the  same  manner  as  those  for  water.  The  continuous  curves  follow  the 
tabulations  for  ice  I at  -7°  by  Schaaf  and  Williams  (1973) (see  labels  and 
legends  on  the  diagram) . The  dashed  lines  connecting  open  triangles 
indicate  the  tabulations  of  Bertie  et  al.  (1969)  for  ice  I at  -173°C, 

This  temperature  is  well  below  the  atmospheric  temperatures  of  cirrus 
cloud  particles  and  snow.  An  extrapolation  with  respect  to  temperature 
is  not  possible  with  present  knowledge.  Furthermore,  atmospheric  ice 
particles  are  mostly  nonspherical  and  irregular  in  shape,  and  a scattering 
theory  for  such  particles  of  finite  size  with  respect  to  the  wavelength 
does  not  exist  except  for  inapplicable  approximations.  Hence  the  attenuation 
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of  far  infrared  and  submillimeter  radiation  through  cirrus  clouds  and 
sncw  can  be  obtained  only  as  a rough  estimate  for  which  an  exact,  knowledge 
of  optical  constants  is  not  essential. 

It  appears  that,  despite  experimental  difficulties,  the  optical  con- 
stants for  liquid  water — at  least  at  room  temperature  or  24°C — are  fairly 
well  determined  in  the  entire  range  up  to  millini^ter  waves.  At  atmospheric 
temperatures  near  0°C,  their  values  may  be  smaller  than  shown  in  Fig.  1 
by  an  unknown  amount  of  the  order  of  a few  percent.  However  because  of 
mutually  compensating  effects  of  such  changes  on  the  scattering  efficiencies 
of  the  smaller  and  larger  particles  in  a polydisperse  aggregate,  the  error 
in  the  overall  extinction  coefficient  may  not  be  significant  if  the  room 
temperature  data  are  used. 

For  convenience  we  show  in  Table  1 the  actual  wavelengths  we  chose 
for  our  calculations  together  with  the  corresponding  complex  index  values 
as  adapted  from  the  original  data  sunmariaed  in  Fig.  1.  A higher  spectral 
resolution  is  not  necessary  in  view  of  the  expected  smooth  variation  of 
polydisperse  scattering  attenuation  with  wavelength.  For  wavelengths  in 
excess  of  X2  mm  we  used  the  Debye  formula  values  a3  in  an  earlier  study 
(Deirmendjian,  1965)  valid  for  G6C  water.  In  the  28  <X<  1000  urn  spectral 
region,  t'.e  optical  constants  shown  are  valid  for  room  temperature  or 
2A°C,  as  obtained  by  experiment,  even  though  cloud  and  raindrops  are 
certainly  at  a lover  temperature.  Since  the  temperature  dependence  of 
these  constants  seems  to  be  uncertain  at  present  we  shall  use  the  24°C 
values  here  without  further  adjustment. 
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3.  Extinction  coefficient  for  polydisper&ions 

The  theory  of  polydisperse  scattering  to  be  used  may  be  outlined  as 
follows.  Let  the  number  density  of  cloud  or  rain  drops  in  a typical 
volume  of  air  be  given  by  a continuous  function  n(r),  where  r is  the  radius 
of  the  spherical  particles.  Then,  if  °ex(m»r)  t^e  extinction  (scattering  + 
absorption)  cross  section  as  a function  or  r and  the  optical  constants 
v and  x,  where  m = v - ix  is  the  complex  index  of  refraction,  it  can  be 


shown  that 


0extx»n»n^r) J “ y*2  <?ex(m,r)n(r)dr 


is  the  ia  lume  extinction  coeffiaientj  8 , for  the  polydi3persion,  where  r 

6X  X 

and  r2  are  the  lower  and  upper  limits  in  the  droplet  radius  r,  and  X is 
the  central  wavelength  of  the  incident  radiation.  This  may  consist  of  a 
narrow  band  of  incoherent  unpolarized  radiation  or  a coherent,  mono- 
chromatic wave  train,  provided  that  the  distance  between  drops  is  large  with 
respect  to  X and  that  the  drops'  space  distribution  is  random  and  nonstationary. 

It  is  well  known  that,  for  spherical  particles  of  finite  size  with 
respect  to  X,  the  behavior  of  a is  highly  dependent  on  r and  m.  In 

the  exact  theory  (see  e.g.  van  de  Hulst,  1957)  one  defines  an  extinction 

2 

efficiency  or  normalized  cross  section  Kfix(ni,x)  = cr  / (nr  ) , given  by  a 
so-called  Mie  series  whose  terms  contain  ratios  of  Ricatti-Bessel  functions 
of  real  or  complex  argument  involving  the  normalized  size  x = 2nr/A  and  m, 

These  series  can  be  readily  evaluated  using  appropriate  programs  for  high- 
speed data  processing  systems  (see  e.g.,  Deirmend j 1 an , et  al,  1961; 

Deirmendj ian , 1969).  Thus  the  extinction  integral  in  (2),  after  a change 
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of  variable, takes  the  form 

( 

3 x 

8 (X,m,n)  = ) /*2  x2n(x)K  (m,x)dx.  (3) 

ex  ' 2 it  / I ex 

From  (3)  it  is  clear  that  for  polydispersions  containing  mainly  finite 
particles  with  respect  to  the  wavelength,  8gx  is  critically  dependent  on 
the  form  and  magnitude  of  n(x) , the  size  and  size-range  of  the  particles, 
as  well  as  on  their  optical  constants. 

Two  limiting  cases  are  wortn  recalling:  the  geometric  optics  or 

large  sphere  approximation  with  X <<  r;  and  the  long  wave  (low  frequency) 

or  small  particle  approximation,  with  r <<  X,  where  the  .chape  of  the 

particle  is  unimportant.  In  the  first  case,  provided  the  spheres  are 

lossy  and  r.  £ 0,  it  is  well  known  by  physical  argument  (van  de  Hulst, 

1957)  and  mathematical  treatment  of  the  Mie  series  (Chylek,  1973)  that 

x co  limit  of  K£x  is  2.  Hence  if  all  the  particles  in  the  poly-’ is  - 

persion  are  in  the  geometric  optics  domain,  we  have  the  well  known  result 

that  8 is  gxven  by 
ex 

r ? 2 

8 [n(r)  ] = 2 rr  f ' r n(r)dr,  X c r;  < + 0 ('♦) 

CX  A 

or  twice  the  sum  of  the  geometrical  cross-sectional  areas  of  all  the 
particles  in  a unit  volume.  Clearly  the  value  of  (1  in  this  case  is 
independent  of  the  wavelength  and  the  optical  constants  but  very  much, 
a function  of  the  density  distribution  n(r).  In  the  nonabsorbing  case, 
where  k m 0,  the  x > oolimit  of  the  scattering  or  extinction  cross 
section  does  not  exist,  as  conjectured  by  Deirmendjian  (1969)  and  rigorously 
demonstrated  by  ChJ'lek  (1973),  and  the  Integration  (9)  is  not  valid. 

However,  since  we  are  dealing  with  the  infrared  spectrum  where  liquid 
water  always  ohows  some  absorption,  this  question  is  not  of  concern  here. 
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In  the  other  limit,  a small  particle  approximation  for  the  extinction 
efficiency  may  be  given  as  a power  aeries  in  the  size  parameter  x (van  de 
Hulst,  1957;  Penndorf,  1962).  In  particular,  as  shown  by  Deirmendjian 
(1963b)  use  of  the  first  term  in  x in  such  a series,  equivalent  to  setting 

Kex(m,x)  - x,  x « 1,  (5) 


Z = (\>2+x2)2  + 4(v2-k2-H), 


may  be  sufficient,  provided  k is  not  small  compared  to  1 and  x 5 O.oi. 
Expression  (5)  represents  also  the  absorption  efficiency  K in  the  Rayleigh 

<1D 

approximation,  in  which  case  K K , K -*-0  for  any  value  of  k no 

clD  6X  SC 

matter  how  small,  provided  x is  sufficiently  small  (Deirmendj Ian,  1969). 


The  extinction  cross  section  itself  may  be  obtained  from  (5)  in  terms  of 
the  volume  v occupied  by  the  particle,  that  is  by  writing 


2 

xr  K =a 
ex  ex 


36nyx_  v 
Z X 


The  extension  of  this  result  to  a polydispersion  follows  easily  by  sub- 
stituting (6)  into  (2)  to  obtain 


where 


3 - 6 , - f V , r <<  A , 

ex  ab  Z A p 2 


P “ 3 


v f 2 r3n(r)c 
*'r. 


is  the  volume  of  all  the  particles  in  the  polydispersion  per  unit  voiume 
of  space.  Note  that  (7)  is  strictly  valid  only  when  the  largest  particle  r. 
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is  itself  much  smaller  than  X,  but  this  restriction  may  be  relaxed  depending 
on  the  nature  of  the  density  function  n(r).  From  (7)  it  follows  that  in 
this  limiting  case  the  form  of  the  density  function  does  not  affect  the 
extinction  if  the  volume  mixing  ratio  remains  constant, 

4.  Drop-size  distribution  functions 

For  the  size  distribution  function  we  shall  use  the  general  form 

n(r)  = a ra  exp  (-brY)  (8) 


first  used  by  this  author  to  model  both  clouds  and  hazes,  as  well  as  rain- 
drops (Deirmendj ian,  1963a, b).  This  function  vanishes  for  very  small  and 
large  radii  and  has  a single  maximum.  The  constants  a,  a,  b,  and  y are 
real  positive  numbers  that  may  be  related  to  measurable  parameters.  For 
example,  for  the  special  case  y - 1,  it  can  be  easily  shown  by  integration 
of  (8)  between  r ■ 0 and<x>  that 


bCi+1 

FS53T  • 1 ' 


(9) 


where  N is  the  total  number  of  particles  per  unit  volume  in  the  distribution. 

Similarly  the  constant  b is  given  by  the  mode  radius  r or  size  of  greatest 

frequency  by  means  of  the  relation  b » a/r  . (We  note  that  Hansen  (1971) 

c 

has  recently  introduced  a distribution  function  which  is  identical  to  (8), 
except  for  a redefinition  or  the  constants  in  terras  of  other  quantities, 
without  additional  physical  justification  in  our  opinion.  The  conversion 
from  one  set  of  constants  to  the  other  is  easily  derived  by  comparing  the 
two  forms.) 


1 Green  and  Lane  (1964)  cite  a 1939  paper  by  S.  Nukiyama  and  Y.  Tanasewa 
that 

who  found/ this  identical  function  represents  the  size  distribution  of 


droplets  produced  by  "small  air  blast  atomizers". 
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Among  the  desirable  properties  of  the  form  (8)  is  the  e vponential 
decay  in  number  density  with  Increasing  drop  size.  This  and  the  nature 
of  the  Mie  functions  result  in  a rapid  convergence  of  the  integrals  of 
type  (3)  in  the  region  defined  by  the  condition 


where  r is  the  radius  corresponding  to  the  upper  limit  x0  of  integration. 
This  behavior  closely  correrponds  to  natural  atmospheric  polydispersione 
where  there  is  a cutoff  of  the  distribution  at  the  upper  end  due  to 
various  processes.  The  function  (8)  also  vanishes  for  very  small  r due 
to  the  power  factor,  again  corresponding  to  natural  conditions.oo  that  the 
Integral  (3)  may  be  replaced  by  a definite  one  between  the  limits  0 andoo 
without  loss  of  applicability. 

In  the  past  our  H,  L,  and  11  haze  models,  and  C.l  cloud  model  for 
fair  weather  cumulus  have  been  found  useful  in  various  applications  of 
radiative  transfer  theory  to  problems  related  to  the  atmospheres  of  the 
earth  and  Venus  (see,  e.g.,  Hansen,  1971;  Thompson  and  Wells,  1971;  Hansen 
and  Hoy  enter,  1974;  Kattawar  and  Plass,  1975).  Obviously  the  mode 
radius  of  4 vim  in  the  C.l  model  is  too  small  and  the  drop  size  range  too 
narrow  to  simulate  a precipitating  type  cloud.  We  shall  now  introduce 
two  new  cloud  density  functions  to  simulate  fog  or  nimbostratus , called 
cloud  C.5,  and  a large  droplet  component,  called  aloud  C.6.  These  choices 
have  been  mainly  suggested  by  measurements,  such  as  delta's  (1961), 
and  others,  taken  within  clouds  that  are  producing  precipitation.  A 
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more  detailed  discussion  and  justification  is  deferred  to  a planned  sepa- 
rate paper, 

The  specific  form  of  our  cloud  C.5  density  function  is 

, , 4 -2r/3  -3  -1  .... 

n(r)  *■  a r e cm  um  , (11) 

with  a mode  radius  of  r of  6 pm.  For  a nominal  density  of  H * 100  droplets 

c 
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per  cm  , this  model  shows  a liquid  water  content  w of  0.297  g m . It  is 

interesting  to  note  that  Hansen's  (1971)  "mean  effective  radius"  of  10.5  pm 

and  "effective  variance"  of  0.143,  corresponding  to  our  C.5  distribution, 

coincide  almost  exactly  with  those  of  Diem's  measured  distributions  for 

cumulus  congestus  and  nimbostratus  clouds  as  reduced  by  Hansen  (op.  ait.). 

For  cloud  C.6  we  adopted  the  formula 

/ - 2 -r  -3  -1  /10. 

n(r;  *=  a r e cm  pm  (12) 

-1  3 

with  r^  - 20  pm.  For  a typical  density  of  10  drop  per  cm  , this  model 

-3 

shows  a water  content  of  0.025  g m . A simple  sum  of  the  density  functions 
(11)  and  (12)  reproduces  fairly  well  two  of  Okita's  (1961)  droplet  spectra 
(for  August  19,  1958;  measured  within  a nimbostratus  cloud.  Their  liquid 
water  contents  are  given  as  0.33  and  0.46  g m respectively  which  compare 

-3 

well  with  our  models'  0.322  g m 

In  general,  the  modeling  of  cloud  drop  spectra  is  relatively  straight- 
forward, since  clouds  usually  are  fairly  stable  for  the  purposes  of  drop- 
let counting  and  sizing  by  modern  instrumentation.  The  phenomenon  of  rain,  on 
the  other  hand,  represents  highly  unstable  conditions  so  that  it  is  diffi- 
cult to  define  a raindrop  density  function  without  specifying  several 
pertinent  variables  and  meteorological  conditions  for  each  case.  The 
easily  recorded  rainfall  rate,  R,  is  an  integral  over  most  such  variables, 
and  its  dependence  on  droplet  growth  and  breakup,  fall  velocities  and 
updrafcs,  distance  below  cloud,  etc.,  is  not  well  known.  Nevertheless, 


any  extinction  of  electromagnetic  flux  by  rain  must  be  interpreted  in 
terms  of  the  "instantaneous"  (really  some  sort  of  time  and  space  average) 
drop-size  distribution  encountered  by  the  beam.  In  radar  meteorology 
these  questions  are  usually  bypassed  by  the  adoption  of  the  so-called  "Z-R" 
empirical  relations  that  seem  to  work  in  operational  situations  (Battan,  1973). 

Here  we  shall  introduce  two  raindrop  density  functions  of  the  type  (8) 
as  examples  of  moderate  and  heavy  rainfall  conditions,  as  suggested  by 
some  recent  drop-size  distribution  and  simultaneous  rainfall  rate 
measurements.  Our  moderate  rainfall  model,  called  rain-l03  is  designed 
to  represent  drop  spectra  just  below  the  cloud  base  when  a precipitation 
rate  of  about  10  ium  hr  1 is  measured  at  the  ground.  The  corresponding 
density  function  is 

n(r)  ■*  a r^  e m ^ mm  \ (13) 

where  r * 1/3  mm.  The  constant  a is  fixed  for  a nominal  N - 1000 
c 
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drops  per  m and  a liquid  water  content  of  0.509  g m . Note  that  the 

model  (13)  distribution  differs  from  the  one  we  introduced  earlier 

(Deii mend j ian , 1969),  called  rain  L (for  light  rain  aloft),  notably  in  the 

choice  of  a bigger  mode  radius,  for  better  agreement  with  Doppler  radar  daLa. 

Our  heavy  rain  model,  called  rain~S03  is  designed  to  represent  raindrop 
spectra  at  ground  level  during  rainfall  rates  of  the  order  of  5C  mm  hr  ^ . 

Its  density  function  is  given  by 


I 


w- 

r 


, . 6 -lOr  -3  -1 

n(r)  a r e m mm 


where  r =0.60  mm, 
c 
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a is  fixed  for  N = 1000  drops  per  m and  w - 2.11  g m 
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The  two  models  are  illustrated  by  the  solid  continuous  curves  in 
Fig,  2,  labeled  "Rain-10"  and  "Rain-50".  For  comparison,  one  of  the 
Doppler  measurements  obtained  by  Caton  (1966)  corresponding  to  a measured 
rate  of  5.6  mm  hr  ^ and  two  of  the  ground-level  measurements  from 
Uataneo  and  Stout  (1968)  are  also  indicated  on  Fig.  1 by  the  labeled 
broken  curves.  Although  the  agreement  is  not  perfect,  our  models  are  seen 
to  be  fairly  representative  of  these  particular  measurements. 

We  emphasize  that  the  above  cloud-  and  raindrop  size  spectral  models 
are  introduced  to  evaluate  the  far  infrared  and  submilliweter  wave  attenuation 
to  be  expected  in  typical  cases.  They  are  not  intended  as  a detailed 
discussion  of  drop  spectra  and  their  relation  to  cloud  microphysics . 

We  do  suggest,  however,  that  linear  combinations  of  these  models,  together 
with  those  of  the  same  type  we  introduced  in  earlier  work,  are.  capable 
of  reproducing  many  of  the  naturally  occuring  water  drop  spectra  and  the 
resulting  scattering  properties. 

5.  The  extinction  coefficients  for  clouds  and  rain 


T'ne  polydisperse  extinction  coefficients  were  computed  on  an  IBM 
370/158  system  for  the  models  described  in  Sec.  4 at  the  wavelengths  and 
optical  constants  of  Table  1 by  the  numerical  techniques  described  pre- 
viously elsewhere  (Deirmendj ian,  1963a,  1969).  The  problem  of  compu- 
tational instabilities  in  the  evaluation  of  Bessel  functions  of  large 
complex  argument  was  overcome  without  difficulty  by  uring  the  double 
precision  arithmetic  and  forward  recursion  techniques  already  described 
(op.  ait.). 
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Selected  values  of  0 and  the  Volume  absorption  coefficient  0 

CX  &D 

(obtained  by  an  integral  similar  to  (3),  are  listed  in  Tables  2 and  3. 

The  difference,  0 - 8 , “ 8_.  is  the  volume  scattering  coefficient t 

as  defined  in  radiative  transfer  formulations.  The  conventional  units 

of  (neper)  km  ^shown,  should  really  be  interpreted  as  a cross-sectional 

7 5 3-1 

area  per  unit  volume,  or  in  this  case  as  cm" (10  cm  ) , such  that  the 

attenuated  flux  F of  a.  plane-parallel  or  collimated  stream  of  electro- 
magnetic energy  is  given  by 


F 


F exp 
o r 


(15) 


where  Fq  is  the  source  energy  and  L is  the  path  length  in  kilometers 
through  a homogeneous  scattering  medium.  In  Tables  2 and  3 the  numbers 
shown  in  parentheses  were  obtained  by  approximations  (4)  and  (7),  as  the 
case  may  be.  The  definite  integrals  involved  in  these  approximations 
may  be  easily  evaluated  from  their  analytic  expressions  in  terms  of 
garnnra  functions,  after  substituting  the  explicit  forms  for  n(r)  in  each 
case  and  integrating  between  the  limits  r - 0 and  oo  . 

For  the  limits  x^(X,r^)  and  of  the  numerical  integration  (3), 

_3 

we  used  the  criterion  that  n(r^)  and  n(r^)  should  be  less  than  10  n(r£) , 

where  n(r£)  is  the  maximum  of  the  density  function.  The  upper 

limit  in  this  case  is  also  governed  by  the  condition  r^  > 2(a+4)/b, 

that  1b  twice  the  value  beyond  which  rapid  convergence  is  assured  according 

to  (10) . The  integration  steps  were  varied  for  each  case  to  assure 

good  accuracy  without  wasting  machine  time.  The  values  shown  in  Tables  2 

and  3 should  be  accurate  to  the  four  significant  figures  Indicated. 
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The  extinction  coefficients  in  Tables  2 and  3 are  displayed  graphi- 
cally in  Fig.  3 where  they  are  plotted  as  solid  dots  logarithmically 
against  the  wavelength.  This  is  marked  in  micrometers  and  millimeters 

on  the  lower  scale  and  the  corresponding  frequency  in  GHz  units  or 
9 -1 

10  sec  is  indicated  on  the  top  scale.  The  ordinate  scale  on  the  left 
shows  the  extinction  coefficient  in  nepers  per  kilometer  and  on  the  right- 
hand  scale  in  decibels  per  kilometer  for  convenience  (conversion  factor 
nepers  to  decibels:  4.343).  The  curves  are  fitted  "by  eye"  to  the 

data  points  and  their  dashed  extensions  connect  the  small  particle 
approximations  shown  in  parentheses  in  Tables  2 and  3.  The  dashed 
straight  lines  with  arrows  pointing  to  the  left  are  the  asymptotic  values 
for  large  spheres  according  to  (4).  For  completeness  and  comparison  we 
also  extended  the  extinction  evaluations  of  our  original  fair  weather 
cumulus  or  cloud  C.l  model  (Deirmendj ian,  1963a,  1969)  to  the  far  in- 
frared with  the  use  of  the  new  optical  constants  shown  in  Table  1.  (The 

_3 

C.l  model  has  a mode  radius  r of  4 pm  and  a water  content  of  0.063  g m 

-3 

for  a drop  density  of  100  cm  .)  The  resultant  infrared  and  submillimeter 
extinction  curve  is  included  in  Fig.  3 fitted  to  the  values  listed  in 
Table  2. 

Tables  2 and  3 list  the  volume  extinction  and  absorption  coefficients 
for  normalized  number  densities  corresponding  to  the  various  models  and 
may  be  used  to  scale  these  coefficients  to  other  number  densities  as 
given  by  observation  or  equivalent  theoretical  models.  It  is  convenient 
and  instructive  to  also  consider  the  extinction  coefficient  per  unit 
maea  of  particles.  This  is  easily  derived  from  the  tables  in  terms  of 
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the  number  of  particles  N per  unit  volume,  per  unit  mass  of  substance. 

‘ m 

From  the  volume  mixing  ratio  in  (7)  and  the  integrated  number  density 
N as  in  (9) , we  have 


N 


m pV 


(16) 


where  p is  the  density  of  the  particle  substance.  In  the  case  of  water 

clouds  and  hazes,  the  liquid  water  content  w - pV  , usually  expressed 

P 

in  grams  per  cubic  meter,  may  be  obtained  from  (16)  by  expressing 

3 3 6 

in  cm  per  m with  p ■ 1 in  which  case  w » 10  V . Consequently  10  N 

p m 

-3 

represents  the  number  of  particles  per  unit  volume  per  g m of  liquid 


water  content  and  the  mass  extinction  coefficient,  y , is  obtained  from 

ex 


C by  wri ting 
ex 


Yex[n(r),w] 


a N 

1°  TT  6 

N ex 


w V 
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(17) 


In  Table  4 we  list  sample  values  of  the  mass  extinction  coeffi- 
cients for  two  of  our  cloud  models,  C.i  and  C,5,  the  rain  model  rain-10, 
and  for  comparison  our  original  haze  model  L (r  “ 0.07  pm)  corresponding 
to  a continental  type  distribution  (Deirmendjian,  1969) • The  use  of 
parentheses  in  Table  4 has  the  same  meaning  as  in  Tables  2 and  3.  The 
reader  can  derive  other  values  by  use  of  equation  (17)  and  the  liquid 
water  contents  for  each  model.  These  values  may  be  used  to  scale  the 
extinction  of  clouds  and  rain  for  liquid  water  contents  other  than  those 
adopted  in  our  models,  provided  that  the  size  distribution  can  be 
approximated  by  these  models.  In  fact,  as  is  clea  from  Table  4,  for  a 
given  water  content  the  amount  of  extinction  produced  by  water  particles 
depends  very  much  on  their  relative  size  nnd  number  distribution. 
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Compare  for  example  the  X ->  0 values,  characteristic  say  of  X0.70  ym 
or  visible  light,  for  rain  and  haze  where  there  is  a thousandfold 
increase  in  extinction  per  unit  mass  of  liquid  water.  The  mass  extinction 
coefficient  at  any  given  wavelength  becomes  independent  of  size  only  when 
all  the  particles  in  the  distribution  are  small  compared  to  the  wavelength 
as  follows  from  equation  (7). 

6.  Discussion  and  ronclusions 

As  may  be  seer  from  Tables  2 and  3 and  Fig.  3,  the  highest  extinctions 

should  be  observed  tn  nimbostratuti-type  clouds  (or  fog)  at  X < 200  ym. 

Here  the  computed  extinction  coefficients  may  even  be  underestimated  by 

factors  of  two  or  three  considering  that  real  clouds  may  contain  200  to 
3 

300  droplets  per  cm  (Mason,  1971).  However,  beyond  X100  ym  the  cloud 
extinction  values  fall  more  or  less  as  X 1 and  are  exceeded  by  those  for 
moderate  rain  at  millimeter  wavelengths. 

A comparison  of  the  extinction  curves  for  the  fair  weather  cumulus 
C,1  and  stratus  cloud  C.5  shows  them  to  be  approximately  parallel,  both 
showing  the  relative  "window"  in  the  10  < X < 12  ym  region.  At  shorter 
wavelengths  toward  the  visible  the  extinction  ratio  8ex(C.l)/8ex(C.5) 
is  given  by  that  of  the  asymptotic  values  for  X -*•  0 listed  in  Table  2, 
that  is,  0.37,  which  is  also  the  ratio  of  the  corresponding  integrated 
geometric  cross  sections  (4).  At  'wavelengths  longer  than  XI  mm  the 
same  ratio  is  seen  to  be  0.21,  which  is  exactly  that  of  the  water  con- 
tents of  the  two  models.  The  contribution  of  cloud  C.6,  as  expected,  is 
seen  to  be  quite  small  compared  to  those  of  the  other  models  except 
in  the  region  X > 100  ym  where  it  is  of  the  same  order  as  the  cloud 


C.l  extinction. 
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In  the  real  world,  of  course,  the  extinction  by  clouds  may  vary 
considerably — though  not  by  an  order  of  magnitude — from  the  estimates 
given  here,  depending  on  the  cloud  type,  water  content,  homogeneity, 
etc.  Our  estimates  are  for  illustrative  purposes  only  and  do  not 
include,  for  example,  the  very  important  contribution  by  the  strong  water 
vapor  absorption  bands  known  to  exist  in  the  entire  region.  Although 
this  subject  is  not  within  the  scope  of  this  paper,  we  note  that, 
according  to  theoretical  estimates  (Birch  et  al,  1969),  the  H^O  absorption 
may  exceed  the  cloud  drop  extinction. 

The  rain  extinction  curves  are  quite  smooth  and  approach  their 
asymptotic  values  quickly  for  wavelengths  shorter  than  100  ym  or  so. 

At  the  other  end  of  the  spectrum,  note  from  the  last  line  in  Table  3 that 
at  A3. 3 cm  the  ratio  7.68  of  the  extinctions  for  the  rain-50  and  rain-10 
models  differs  significantly  from  that  of  their  respective  water  con- 
tents, namely  2.11/0.509  = 4.15.  This  demonstrates  that  the  small  sphere 
approximation  (7)  is  not  good  enough  at  this  wavelength,  as  pointed  out 
elsewhere  (Delrmendjian,  1963b).  The  maximum  in  rain  extinction  occurs 
somewhere  between  A2  and  A5  ran  for  both  models,  in  a region  where  the 
cloud  effect  is  small  and  dropping  rapidly  toward  longer  waves. 

A rough  comparison  of  our  rain  extinction  estimates  with  indepen- 
dent evaluations  may  be  obtained  by  graphical  interpolation  on  data 
adapted  from  the  literature  by  Dudzinsky  (1974).  For  a rainfall  rate  of 
50  mm  hr  \ our  value  and  that  estimated  by  D.  E.  Setzer  (Dudzinsky,  1974, 

Fig.  1)  coincide  exactly  at  A10  mm,  while  at  A33  mm  our  value  is  somewhat 
lower  than  Setzer' s estimate.  Similarly,  our  rain-50  model  values  fall 
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between  the  "world"  measured  maximum  and  minimum  attenuations  as  given  by 
R.  G.  Medhurst  and  summarized  by  Dudzinaky  (1974,  Fig.  3).  For  wave- 
lengths shorter  than  X5  mm  there  seem  to  be  few,  if  any,  observations 
and  our  present  theoretical  values  may  be  useful  in  the  design  of 
future  experimental  programs  and  applications. 

For  the  cloud  extinction  values  in  the  submillimeter  region  also 
there  are  few  Independent  estimates  and  experimental  values  to  allow  a 
comparison.  A good  determination,  on  the  basis  of  a series  expansion  of 
Kex>  is  provided  by  Platt  (1970).  Platt  made  use  of  the  first  three 
terms  of  the  power  series  in  x,  instead  of  only  the  first,  as  in  equation 
(5) , to  obtain  the  extinctior  contribution  from  the  larger  droplets  at 
X200  and  1337  um,  his  two  shortest  wavelengths,  together  with  the  best 
optical  constants  available  to  him.  We  compare  below  his  determinations 
for  a typical  fog  distribution  with  our  own  model  values  from  Table  4, 
expressed  in  units  of  nepers  par  kilometer  per  unit  liquid  water  content: 


X 200  um 

X 337 

Platt  (1970,  Table  5) 

21.4 

9.46 

Cloud  C.l  (N  - 1599) 
m 

15.6 

9.29 

Cloud  C . 5 (N  - 337) 
m 

18.9 

9.93 

Considering  the  differences  in  method  and  drop-si  jtribution,  the 
agreement  is  quite  good,  especially  at  the  longer  wavelength  as  ex- 
pected, and  corroborates  Platt's  (1970)  method  of  approximation  in  these 


cases . 
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Finally  a comparison  with  available  measurements  of  attenuation  in 
the  natural  atmosphere  should  be  made  as  a check  of  the  validity  of  our 
models.  Unfortunately  there  are  too  few  published  measurements,  over 
a wide  enough  portion  of  the  far  infrared  spectrum,  for  an  adequate 
comparison.  The  point  marked  "Fog  (obs)"  on  Fig.  3 is  adapted  from  a 
determination  at  the  HCN  maser  wavelength  of  A337  urn  (Burroughs  et  al, 

1966)  obtained  in  fog  with  a visibility  of  35  m.  This  is  equivalent  to  an 
extinction  coefficient  of  about  96  nepers  per  kilometer  at  visible 
wavelengths  according  to  a modification  of  Koschmieder 's  formula 
(Middleton,  1952),  or  about  double  the  short  wave  asymptotic  value  for  our 

cloud  C.5  model.  This  implies  that  the  fog  contained  some  200  droplets 

3 -3 

per  ci  and  a liquid  water  content  of  about  0.6  g m , not  at  all  unrea- 
sonable numbers  according  to  existing  measurements  (Mason,  1971).  Our 

corresponding  extinction  coefficient  at  A337  pm  from  Table  2 or  4 is 
_ 1 

5.9  km  which  checks  very  well  with  the  Burroughs  et  al,  (1966)  measurement 
of  5.8  km  1 plotted  on  Fig.  3 and  constitutes  a kind  of  validation  of  our 
model  to  represent  some  fogs. 

The  point  on  Fig.  3 marked  "water  vapor  (obs)"  is  adapted 
from  the  same  source.  It  represents,  a water  vapor  absorption  equivalent 
to  an  attenuation  of  11.5  km  \ measured  over  a horizontal  path,  indicating 
the  magnitude  of  this  component  in  clear  air  in  the  absence  of  fog  and  with 

-3 

an  absolute  humidity  of  4.84  g m . Within  precipitating  clouds  or  fog 
under  saturation  conditions,  the  absolute  humidity  may  be  of  the  order 

-3 

of  10  g m with  a corresponding  increase  in  water  vapor  absorption- 
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An  idea  of  the  magnitude  of  attenuation  to  te  expected  under  ex- 
treme conditions  of  both  fog  and  heavy  rain  in  the  path  of  a beam  of 
X337  pm  radiation  is  given  by  the  crossed  circle  on  Fig.  3,  marked 
’’total  estlm.  arrived  at  as  follows: 


Nepers  km  ^ 

Decibels  : 

H^O  absorption  (10  g m-3) 

24.0 

104. 

-.3 

Cloud  C.5  (w  » 0,6  g m ) 

5.90 

25.6 

Cloud  C.6 

0.68 

2.95 

Rain-50 

4.06 

17,6 

Total  34.64  150.2 


In  this  example  we  see  that  rain  accounts  for  only  about  12  percent  of  the  total, 

the  heaviest  contribution  being  that  of  water  vapor  absorption. 

As  to  attenuation  by  other  types  of  cloud  and  precipitation  com- 
posed of  ice  particles,  snow,  and  hall,  ar  accurate  theoretical  estimate 
is  unobtainable  at  present,  since  (a)  the  effects  of  the  irregularity 
of  their  shape  are  not  well  known  except  for  highly  idealized  models, 
and  (b)  the  optical  constants  for  ice  in  the  submillimeter  region  have 
so  far  been  determined  only  for  unrealistically  low  temperatures. 

However,  a rough  idea  of  the  attenuating  properties  of  such  particles 
could  b et  obtained  in  certain  cases  by  using  the  applicable  asymptotic 
approximations  discussed  earlier. 

The  scope  of  the  present  study  should  be  sufficient  to  provide 
good  estimates  of  the  attenuating  effects  of  terrestrial  water  clouds 
and  rain  for  wavelengths  up  to  X20  mm.  Our  results  apply  only  to 
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power  degradation  due  to  incoherent  scattering  but  not  to  phase  and  co- 
herence losses  caused  by  other  types  oC  scattering  and  turbulence.  Water 
vapor  absorption  under  atmospheric  conditions  needs  further  study.  There 
seem  to  be  several  window-like  bands  in  the  region  0.83  < X < 8.6  mm, 
according  to  papers  presented  in  two  recent  symposia  (Fox,  1971;  Button, 
1974)  based  on  theoretical  estimates  and  laboratory  determinations. 

From  our  present  results  it  would  thus  seem  that,  in  the  presence  of 
nonprecipitating  water  clouds,  the  presence  of  a relative  "window"  around 
XI. 3 mm  or  a frequency  of  230  GHz  might  be  postulated.  At  this  frequency 
Wrixon  (1974)  obtained  a fairly  good  estimate  of  the  clear  sky"  zenith 
attenuation,  of  about  1 dR  or  an  optical  thickness  of  0.23,  by  using  the 
reliable  Bouguer-Langley  method  and  solar  millimeter  radiation.  This 
is  of  the  same  magnitude  as  that  of  1 km  of  cloud  C.l.  Precipitating 
type  clouds,  fog,  or  rain,  however,  would  increase  the  value  by  large 
factors  depending  on  their  th'c  .ness  and  direction  of  "line  of  sight". 

Wg,  of  course,  have  also  generated  the  angular  scattered  intensity 
and  polarization  parameters  corresponding  to  the  new  cloud  and  raindrop 
density  models  introduced  here.  Their  analysis  and  implications  in  cloud 
mlcrophysice  research  will  be  discussed  in  separate  papers. 
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Table  1 

WAVELENCTH  AND  COMPLEX  INDICES  OF  REFRACTION 


A 

m 

X 

m 

12  .pm 

1.111 

- 

0.1991 

500. pm 

2.22 

0.7401 

17.  " 

1.376 

- 

0.4291 

700.  " 

2.32 

- 

0.8901 

28.  " 

1.549 

- 

0.3331 

1000,  " 

2.50 

- 

1.091 

40.  " 

1.519 

- 

0.3851 

60.  " 

1.703 

- 

0.5871 

2 .mm 

2.5604 

- 

0.89471 

100.  " 

2.06 

- 

0.5511 

5.  " 

3.1918 

- 

1.76571 

140.  " 

2.07 

- 

0.4701 

10.  " 

4.2214 

- 

2.52591 

200.  " 

2.08 

- 

0.5091 

20.  " 

5.8368 

- 

3.00461 

337  . " 

2.20 

- 

0.6001 

33.  " 

7.1755 

2 . 86421 
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Values  from  an  earlier  run  with  m*l. 369-0. 4381 
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Table  3 


RAIN  VOLUME  EXTINCTION  AND  ABSORPTION  COEFFICIENTS 
(Neper  kin  \ N**10Bm 


X 

Rain 

-10 

Rain 

-50 

^ex 

Bab 

!ex 

Bab 

(X-O) 
60  .pm 

(1.310) 

1.396 

0.598 

(3.519) 

100.  " 
140.  " 

1.434 

1.462 

0.607 

0.635 

3.770 

1.566 

o 

o 

c* 

1.502 

0,666 

3.905 

1.677 

337  . " 

1.577 

0.715 

4.062 

1.766 

500.  " 

1.651 

0.769 

30 

c <1 

1.868 

700.  " 

1.726 

0.816 

4.383 

1.957 

1000.  " 

1.823 

0.922 

4.555 

2.188 

2 .mm 

1.950 

1.017 

5.043 

2.518 

5." 

1.273 

0.753 

5.015 

2.548 

10,  " 

0.3391 

0.2640 

2.293 

1.452 

20.  " 

0.0675 

0.0632 

0.5101 

0.4470 

33.  " 

0.0193 

0.0188 

0.1482 

0.1408 
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Table  4 


MASS  EXTINCTION  COEFFICIENTS  FOR  HAZE, 

-l  -3 

(Yex  in  Neper  km  ™ per  gm  liquid 

CLOUDS,  AND  RAIN 
water  content) 

X 

Haze  L 

Cloud  C.l 

Cloud  C.5 

Rain-10 

w**1.167'10~3»m"3 

w=0 .06255  gm"3 

„ “3 

w=0.2969  Rm 

w=*0.5091  gm 

(A-*0) 

(3117.) 

(250.1) 

(142.8) 

(2.573) 

16.6pm 

247.6 

17.0  " 

257.8 

168.3 

100 . " 

(36.8) 

43.01 

69.72 

2.816 

200.  " 

(16.8) 

15.56 

18.92 

2.950 

337  . " 

(10.5) 

9.293 

9.932 

3.097 

500.  " 

( 7.21) 

7.301 

7.527 

3.243 

1 .mm 

( 3.87) 

(3.87) 

3.924 

3.580 

2.  " 

( 1.60) 

(1.60) 

(1.60) 

3.830 

32. 
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